The problem of a metal impurity in an oxide: ab initio study of electronic and 
structural properties of Cd in Rutile Ti0 2 



o 
o 

Oh! 
< 
00 



o 

CO 

cd 



T3 

o 

(N 
> 

o\ 
o 

o 
m 
o 

-I— > 
cd 



i 

c 

o 
o 



X 
S3 



Leonardo A. ErricoQ Gabriel Fabricius^, and Mario RenterfgQ 

Departamento de Fisica, Facultad de Ciencias Exactas, 
Universidad Nacional de La Plata, C.C. N° 67, 1900 La Plata, Argentina 
(Dated: February 2, 2008) 

In this work we undertake the problem of a transition metal impurity in an oxide. We present an 
ab initio study of the relaxations introduced in TiC>2 when a Cd impurity replaces substitutional^ 
a Ti atom. Using the Full-Potential Linearized- Augmented-Plane- Wave method we obtain relaxed 
structures for different charge states of the impurity and computed the electric-field gradients (EFGs) 
at the Cd site. We find that EFGs, and also relaxations, are dependent on the charge state of the 
impurity. This dependence is very remarkable in the case of the EFG and is explained analyzing the 
electronic structure of the studied system. We predict fairly anisotropic relaxations for the nearest 
oxygen neighbors of the Cd impurity. The experimental confirmation of this prediction and a brief 
report of these calculations have recently been presented [P.R.L. 89, 55503 (2002)]. Our results for 
relaxations and EFGs are in clear contradiction with previous studies of this system that assumed 
isotropic relaxations and point out that no simple model is viable to describe relaxations and the 
EFG at Cd in Ti02 even approximately. 

PACS numbers: 71.55.Ht,61.72.-y,71.20.Ps 



I. INTRODUCTION 



The problem of metal impurities in oxides is a chal- 
lenge in solid-state physics both from a fundamental and 
applied point of view. Description of impurity levels in 
oxide semiconductors has a great technological interest, 
but a complete theoretical description of the problem is 
very difficult. In effect, impurities introduce atomic re- 
laxations in the host and modify the electronic structure 
of the system, being the interplay between these two ef- 
fects one of the principal difficulties in the theoretical 
approach. Some experimental techniques used in mate- 
rial science introduce probe-atoms into the systems to 
be studied?^ giving valuable physical information seen 
by the probe-atom, usually an impurity in the system. 
In particular, TDPAC spectroscopy 3 is a hyperfine tech- 
nique that enables a high resolution determination of 
the EFG tensor at impurity sites. Due to the r~ 3 de- 
pendence of the EFG operator from the charge sources, 
the EFG is mostly originated in the electronic charge 
density close to the impurity nucleus, which in turn re- 
flects the probe chemistry with its neighborhood. Thus, 
it would be important that the theoretical approach to 
the metal-impurity problem would be able to predict 
the EFG tensor with the required accuracy. Since the 
EFG is very sensitive to the anisotropic charge distri- 
bution close to the nucleus, for its accurate calculation 
the entire electronic configuration of the host, perturbed 
by the presence of the impurity, has to be determined. 
This kind of calculations should be performed at the ab 
initio level and can be done in the frame of the Den- 
sity Functional Theory (DFT). A well proven method to 
solve the all-electron DFT equations in solids is the Full- 
Potential Linearized- Augmented-Plane- Wave (FLAPW) 
method. - This method was able to predict with high ac- 
curacy and precision the EFG at undisturbed lattice sites 



in a large variety of pure systems However, on sys- 
tems with impurities, very few calculations have been 
performed and the state-of-the-art is far from being rou- 
tinely in this &e\dJ^^JJ^L1^^5 Tne first fully self- 
consistent ab initio determination of the EFG tensor at 
an impurity site (Cd) in an oxide (Ti02) has recently 
been reported^. In that work we performed a FLAPW 
calculation of the relaxations introduced by the impurity 
and studied their interplay with the electronic structure 
of the system, predicting highly anisotropic relaxations 
of the nearest neighbors of the impurity and a drastic 
change in the orientation of the principal component of 
the EFG tensor. This prediction was confirmed in the 
same work by a key PAC experiment. Due to the com- 
plexity of the physical situation we want to describe, even 
in the framework of ab initio calculations some variables 
have to be determined through comparison with exper- 
iments. For example, we performed the calculations for 
two charge states of the impurity and one of them could 
be discarded because it gave a value for the EFG asym- 
metry parameter, 77, incompatible with experience. 

The central purpose of this work is to present and dis- 
cuss the electronic structure of the Ti02:Cd impurity sys- 
tem with more detail that we could do in Ref. Here, 
we study the relationship between the calculated EFG at 
the impurity site - the key experimental quantity - and 
the interplaying physical quantities: the structural relax- 
ations and the character and filling of the impurity state. 
We also discuss in detail the precision of our calculations 
studying its dependence on the different parameters that 
control the precision (impurity dilution, base dimension, 
average in k-space, etc.) and its dependence on the ap- 
proximation used for the exchange-correlation potential. 

Another important question to be addressed concern- 
ing all experimental techniques that introduce impurity- 
tracers in solids is at what extent simple models usually 
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used to predict the measured quantities are acceptable 
in view of the modifications introduced by the presence 
of the impurity. Concerning the EFG, simple approxi- 
mations like the point-charge model 2 with antishielding 
factor o 17 i 18 and isotropic relaxations of nearest-neighbor 
atom o 18 i 19 are commonly used. In this work we give ar- 
guments against the applicability of these simple models 
for the studied system. 

This paper is organized as follows. In section II we de- 
scribe the method of calculation used in the present work 
and include some preliminary calculations for an unre- 
laxed cell in order to discuss methodologically the point 
of the charge state of the impurity. In section III.(A-C) 
we present and discuss the results of performing relax- 
ations of the Cd nearest-neighbors (NN) in a 72-atoms 
super-cell for two different states of charge of the impu- 
rity, while in section III.D we study the accuracy of these 
results performing several additional calculations. Even 
if section III.D could be found rather technical by people 
not specially interested in first-principles calculations, we 
think it will be useful for the ab initio community inter- 
ested in the problem of impurities in solids. In section 
III.E we compare our results with previous studies of this 
system. Finally, in section IV we present our conclusions. 



II. METHOD OF CALCULATION 

A. A general scheme 

Our aim is to calculate from first principles the struc- 
tural relaxations produced in the Ti02 lattice when a Cd 
impurity replaces a Ti atom and the electronic structure 
of the resulting system. In particular we are interested 
in the EFG tensor at the Cd site. The general scheme 
we adopted to deal with this problem was the follow- 
ing: we considered a super-cell (SC), containing a single 
Cd impurity, repeated periodically and performed first- 
principles calculations We studied the relaxation intro- 
duced by the impurity computing the forces on the Cd 
neighbors and moving them until the forces vanished. 
The calculations were performed with the wien97 im- 
plementation, developed by Blaha et al.,~ of the Full- 
Potential Linearized- Augmented-Plane- Wave (FLAPW) 
method and we worked in the LDA approximation using 
the Perdew and Wang parametrization for the exchange- 
correlation potentials 21 In this method the unit cell is 
divided into non-overlapping spheres with radii Ri and 
an interstitial region. The atomic spheres radii used for 
Cd, Ti and O were 1.11, 0.95 and 0.85 A, respectively. 
We took for the parameter KKmax, which controls the 
size of the basis-set in these calculations, the value of 
6 that gives basis-set consisting in 4500 LAPW func- 
tions for the SC described in section Hi Bl We also intro- 
duced local orbitals (LO) to include Ti-3s and 3p, 0-2s 
and Cd-4p orbitals. Integration in reciprocal space was 
performed using the tetrahedron method taking an ade- 
quate number of fc-points in the first Brillouin zone. Once 



self-consistency of the potential was achieved, quantum- 
mechanical-derived forces were obtained according to Yu 
et aZiSS, the ions were displaced according to a New- 
ton damped scheme 2 ^ and new positions for Cd neigh- 
bors were obtained. The procedure is repeated until the 
forces on the ions are below a tolerance value taken as 
0.025 eV/A. At the relaxed structure, the elements 
of the EFG tensor are obtained directly from the Vzm 
components of the lattice harmonic expansion of the self- 
consistent potential (a more detailed description of the 
formalism used to compute the EFG tensor may be found 
in Ref.i). 

There is still an important point to take into account 
concerning the calculation of the electronic structure of 
an impurity system: the charge state of the impurity. 
We let this point to the end of this section as it will 
be easily handled after the analysis of some preliminary 
calculations. 



B. Cell and super-cell 



(a) z 




FIG. 1: (a) Unit cell of rutile T1O2. All the results discussed 
in this article are referred to the axis system indicated in this 
figure assuming that Cd replaces the black Ti atom. Even 
when Ol and 02 atoms are equivalent in Ti02, it will not be 
the case in the impurity system; (b) 72-atoms super-cell used 
in our calculations. Some O and Ti atoms are not shown for 
clearness. 
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The rutile (TiC^) structure is tetragonal (a = b = 
4.5845i A, c = 2.9533i A). The unit cell (shown in Fig. [J) 
contains 2 metal atoms (Ti) at positions 2a (0, 0, 0) and 
(1/2, 1/2, 1/2) and 4 anions (O) at positions 4/ ±(u, u, 
0; u+1/2, 1/2-u, 1/2), with u = 0.30493 7 ^ The super- 
cell considered in the present work consists of twelve unit 
cells of TiC>2 where one Ti atom has been replaced by 
a Cd atom. The resulting 72-atom SC (called 72A-SC 
in the future) has dimensions a' =2a, b' =26, c' —3c 
and is also tetragonal with d /a! =0.97 giving an almost 
cubic lattice. This SC keeps Cd atoms as far as possible 
from each other for the given cell volume. For checking 
purposes we have also considered SCs containing eight 
and sixteen unit cells. We assume that relaxations to be 
performed preserve point group symmetry of the cell in 
this initial configuration. Symmetry restricts 01 and 02 
displacements to yz plane and x direction respectively. 
In order to check this assumption and the stability of the 
obtained solution, at the end of the relaxation process 
we performed new calculations with 01 and 02 atoms 
displaced from their symmetry positions and verified that 
these solutions are a minimum for each system studied. 

The self-consistent calculations were performed taking 
8 fc-points in the irreducible Brillouin zone for the metal- 
lic system (situation (i)), and 2 fc-points for the non- 
metallic ones (situation (ii)). In order to plot the DOS 
we calculate eigenvalues at a denser mesh of 36 fc-points. 



C. Preliminary study: Charge state of the impurity 

In order to discuss the problem of the charge state 
of the impurity we have done some preliminary calcu- 
lations. We have calculated the self-consistent electronic 
structure of the 72A-SC with all atoms in their initial un- 
relaxed positions. We want to analyze the changes in the 
electronic structure of the system caused by the presence 
of the impurity, neglecting for the moment the problem 
of structural relaxations. The density of states of this 
system is compared with the one of pure Ti02 in Fig. [5] 
Pure Ti +4 02~ 2 is a wide band-gap semiconductor with 
the O-p band filled and the Ti-d band empty. Since Cd 
valence is 2+, when a Cd atom replaces a Ti atom in the 
SC the resulting system is metallic because of the lack 
of two electrons necessary to fill up the O-p band. The 
question that arises at this point is if the real system we 
want to describe provides the lacking two electrons (via 
an oxygen vacancy, for example) or not, and if this point 
is relevant for the present calculation. We will show first 
that this point is absolutely relevant for the present cal- 
culation. Comparison of Fig. [UJa) and (b) shows that 
the presence of the Cd impurity in the SC introduces 
the appearance of Cd-d levels and impurity states at the 
top and the bottom of the valence band in the corre- 
sponding DOS. The ones at the top of the valence band 
can be better seen looking at the band structure in this 
energy range shown in Fig. [3J The two bands that are 
immediately above and below the Fermi level when cross- 
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FIG. 2: (a) DOS for pure Ti0 2 . (b) DOS for the 72-atom 
super-cell (unrelaxed). Energies are referred to the Fermi 
level. Note that the Cd-d band is described as very sharp 
peaks indicating that interaction between Cd atoms of dif- 
ferent cells is quite small for the considered super-cell. The 
arrows indicate impurity states in the valence band. 



ing A-point correspond to impurity anti-bonding states 
that are spatially located at Cd and at their 01 and 02 
nearest-neighbor atoms. In particular, the wave function 
of the impurity state that remains almost completely un- 
occupied has character Cd-d yz , 01-p y , and 01-p z . Then, 
providing two electrons to the system implies a drastic 
change in the symmetry of the electronic charge distri- 
bution in the neighborhood of the impurity. Therefore, 
different nearest-neighbors relaxations and also different 
EFG could be expected for different charge states of the 
impurity. 

In the present work, we present calculations for two 
charge states of the impurity system, corresponding to 
two different physical situations: (i) Cd (neutral im- 
purity state), corresponding, e.g., to an extremely pure 
crystal at low temperatures, and (ii) Cd~ (charged ac- 
ceptor state), the system provides two electrons via an 
oxygen vacancy, donor defects, etc. To study case (i) we 
used the 72A-SC already described, to study case (ii) we 
also used the 72A-SC but performed self-consistent cal- 
culations adding two electrons to the SC that we compen- 
sate with an homogeneous positive background in order 
to have a neutral cell to compute total energy and forces 
(this procedure is implemented since version wien97.9 of 
the FLAPW package). We have also simulated situation 
(ii) with an alternative procedure: we replaced in the 
72A-SC the two most distant oxygen atoms by two fluo- 
rine atoms. In this way we provide two electrons to fill 
up the O-p band without introducing any artificial back- 
ground. We expected that the difference between flourine 
and oxygen potentials should modify only slightly the re- 
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FIG. 3: Band structure for the 72-atom super-cell (unre- 
laxed) . The zero of energy is defined as the average Coulomb 
potential in the interstitial region. 



suits. In summary, we performed self-consistent FLAPW 
calculations for the following systems: 

(t) Cd°, (Ti0 2 ) 23 Cd0 2 
(it) Cd-(2e), (Ti0 2 ) 23 Cd0 2 + 2e~ 

Cd~ (fluorine), (Ti0 2 ) 23 CdF 2 



III. RESULTS AND DISCUSSION 



A. Structural relaxations 



TABLE I: Final coordinates of the Cd nearest oxygen neigh- 
bors for the different calculations performed with the 72A-SC 
compared with the ones of pure Ti02- d(Cd-Ol) and d(Cd- 
02) are the distances (in A) from Cd to Ol and 02 atoms, 
respectively. 8i is the angle (in degrees) between z axis and 
foi - r C d- 





d(Cd-Ol) 


d(Cd-02) 


6i 


Ti0 2 


1.944 


1.977 


40.47 


Cd° 


2.153 


2.108 


39.59 


Cd"(2e) 


2.185 


2.111 


39.55 


Cd~ (fluorine) 


2.191 


2.121 


39.73 



with this problem. 

Anisotropy in the relaxations of the nearest oxygen 
neighbors of the Cd impurity can be understood by in- 
spection of Fig. Q] Stretching of Cd-02 bond implies 





Ji Cd 



FIG. 4: Planes XZ and YZ (see Fig.rjJ containing 02 and Ol 
atoms, respectively, with their neighbors. The arrows indicate 
the displacement of the oxygen atoms from the unrelaxed to 
the final relaxed positions in the 72A-SC[Cd~ (2e)] system. 
The size of the relaxation has been duplicated in order to 
better visualize the effect. 

a considerable shortening in Ti-02 bonds. However, 
stretching of Cd-Ol bond doesn't affect so much Ti-Ol 
bonds since the structure is more open in this direction. 
So, at the end of the relaxation, Cd-Ol bond stretches 
almost twice than Cd-02 bond. 



Let us first consider the relaxation of only the 6 near- 
est oxygen neighbors of the Cd impurity (Ol and 02 in 
Fig [TJ) In Table U we compared the results of the re- 
laxation of these oxygen atoms for the different systems 
studied. We see that for both charge-states of the impu- 
rity the relaxations are quite anisotropic, with the Cd-Ol 
distance larger than Cd-02 distance, opposite to the ini- 
tial unrelaxed structure. This result is opposite to what 
other authors have assumed in previous studies of this 
syste m 18 ' 19 and confirms the tendency predicted in our 
previous calculation with a much smaller SC. 25 As it can 
be seen in Table U the difference in the charge state of the 
impurity affects essentially the relaxation of 01 atoms 
that present a slightly larger relaxation for the charged 
impurity. Relaxation for the two ways of simulating the 
charged state of the impurity (case Cd - ) are very similar 
indicating that both approaches are well suited to deal 



B. Electronic structure 

In Fig. [5] we show the bands for the 72A-SC for the two 
states of charge considered for the impurity. Comparison 
of Fig. OJa) with Fig. [3] shows that as a consequence of 
relaxation the outermost impurity level falls down in en- 
ergy and goes into the O-p band becoming half-occupied. 
Relaxation increases the Cd-Ol distance and this pro- 
duces a softening of the Cd-d— Ol-p interaction and the 
fall down in energy of the anti-bonding impurity states. 
Figure [5] shows that the band structure of the neutral and 
charged relaxed structures are very similar, but the out- 
ermost impurity state is a slightly raised in energy when 
it is completely filled (case (ii)) as a consequence of the 
larger Coulomb repulsion. Comparison of Fig. [Sfa) and 
Fig. El^c) shows the little fall down in energy mentioned 
for the anti-bonding impurity states an also a little rise 
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FIG. 5: Band structure of the relaxed 72A-SC for different charge states of the impurity: (a) Cd (neutral state) and (b) 
Cd~(2e) (charged state). The zero of energy as in Fig. 3. 



for the bonding ones at the bottom of the O-p band. A 
shift upwards of about 1.5 eV of Cd-d levels from the 
un-relaxed 72A-SC(0) to the relaxed one is also present. 
To look at the orbital composition of the impurity states 
we plot in Fig. [Tfa) and (b) the PDOS for Cd-d, 01- 
p and 02-p for the two charge states of the impurity. 
The impurity state near the top of the valence band has 
Cd-d yz , Ol-py and 01-|> z character and, as we men- 
tion looking at the bands, it is shifted upwards in energy 
when it is completely filled (case Cd - ). In Fig. [7|a)(Cd 
atom) and (b)(Cd atom) it can also be seen the pres- 
ence of other impurity state with components Cd-t4.2_.y2 
and Cd-c?3 Z 2_ r 2 which is completely filled in both cases. 
This impurity state is located mainly at Cd-d x 2 _ y ? and 
02-p x but involves also contributions from 01-p z , Cd- 
d 3z 2_ r 2 and 03-p x (03 is the NN of 02 atom in the x 
direction). From the present calculations the outermost 
impurity state has an occupation of around 1.3 e for the 
72A-SC(Cd°) and 2 e for the 72A-SC(Cd-). We identify 
this 0.7 additional electrons in this state as the driving 
force that produce the slightly larger relaxation for 01 
atoms in the charged cell with respect to the neutral cell. 

We want to mention that the fact that the impurity 
level falls at the Fermi energy in the Cd° case is not for- 
tuitous. Due to Coulomb repulsion, the impurity level 
would fall below the Fermi energy if it was empty an 
above if it was filled being the only self-consistent solu- 
tion to be half-occupied at exactly the Fermi energy. The 
occupation of the impurity level in the 72A-SC(Cd°) is 
therefore a constant number and fairly independent of 
small fluctuations of charge in the cell. 



C. Electric-Field Gradients 

In Table [Til w e show the results for the Vu principal 
components of the EFG tensor for the three systems stud- 
ied. The resulting EFGs for the two approaches used to 
simulate the charged impurity are very similar, the dif- 
ference of 0.4 x 10 21 V/m 2 in components Vxx and Vyy 
is within one could expect for the small difference found 
in the oxygen positions, since EFG is very sensitive to 
small structural changes. These results agree very well 
(in magnitude, symmetry and orientation, see Table [TTJ 
with the experimental results obtained for the EFG at 
Cd impurities substitutional^/ located at cationic sites in 
rutile Ti02-— The difference between the EFGs obtained 
for the charged and neutral cells is very remarkable: the 
sign, direction and absolute value of the largest Vu com- 
ponent (V33) are different in both situations, and also the 
value of the asymmetry parameter rj. The high rj value 
obtained for 72A-SC(0) shows that the electron availabil- 
ity present in the sample leads the impurity to be in a 
charged state. 

In order to investigate the origin of the difference in the 
EFG for the two charge states of the impurity we con- 
centrate in the valence contribution to the EFG which 
originates in the asymmetry of the valence charge distri- 
bution inside the muffin-tin sphere. The valence con- 
tribution is usually dominant in FLAPW calculations 
and can be split in the different orbital symmetries.- In 
Table IHII we show the total valence contribution to Vu 
and its components arising from p and d orbital symme- 
tries. We see that the largest differences correspond to 
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TABLE II: EFG tensor principal components at Cd site, Vn(in 10 21 V/m 2 ), for the relaxed structures of the different systems 
considered in our calculation compared with experiments and the calculation of Sato et al. r\ — (Vu — V22)/V33 ( | V33 1 > | V22 1 > 
|Vii|). In the last row EFG tensor refers to Ti site in pure Ti02- Q=0.83 b(Q=0.24 b) was used to calculate V33 from the 
experimental quadrupole coupling constant vq at 111 Cd( 49 Ti) sites. 



Vxx Vyy Vzz V33 V33-direction r) 

72A-SC(0) -7.16 +6.82 +0.34 -7.16 X 0.91 

72A-SC(-2) -2.87 +4.55 -1.68 +4.55 Y 0.26 

72A-SC(F) -2.46 +4.10 -1.63 +4.10 Y 0.20 

Exp.iS 5.23(5) ... 0.18(1) 

Exp.^ 5.34(1) ... 0.18(1) 

Exp. (singlecrystal)i£ 5.34(1) X or Y 0.18(1) 

C&\cM +1.54 +3.56 -5.09 -5.09 Z 0.39 

Exp. (pure Ti0 2 )^ 2.2(1) Z 0.19(1) 
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FIG. 6: DOS of the (a) relaxed 72A-SC (Cd°), (b) relaxed 
72A-SC[Cd-(2e)], and (c) unrelaxed 72A-SC (Cd°). Figure 
(c) is the same as Fig. [2fb) and has been repeated here in 
other scale for the sake of comparison. Energies are referred 
to the Fermi level. 



d-components of Vu- This difference is originated in the 
filling of the impurity state at the Fermi level that has 
an important component of Gd-d yz symmetry as can be 
seen in Figs.[7£a) and (b)(Cd atom). A simple analysis in 
terms of partial charges 2 ^ shows that the effect of adding 
Sn electrons to an orbital d yz is to produce a change in 
Vu components given by: SVxx= Id^n, 5Vyy= - Id^ri/2, 
SVzz— - where Id is proportional to (1/r 3 ) for d 

orbitals inside the muffin tin sphere, fntegration of unoc- 



cupied C&-dy Z PDOS from Fig. [7ja) (Cd atom) gives 5n= 
0.074. Inspection of Table IIIII shows that the changes 
SVu are quite well described by this estimation giving for 
Id a value around 83 x 10 21 V/m 2 . Another interesting 
point is the presence of d contributions to Vu in the 72A- 
SC[Cd~(2e)] although the Cd-d band is completely filled 
(see Table IIIII) . This curious point may be explained be- 
cause radial d wave functions are energy dependent, so, 
electrons in anti-bonding impurity states (at the Fermi 
level) contribute to Vu with a different (l/r 3 ) factor than 
the ones of the same symmetry in Cd-d levels (between 
6 and 7 eV below the Fermi level). We have found an in- 
crement of 30% in (l/r 3 ) when going from an energy of 
-7 eV to the Fermi level. This dependence of (l/r 3 ) with 
energy was also mentioned and verified by Blaha et al. in 
Ref. [29| for Cu in CU2O. Finally, it has to be mentioned 
that the difference in p contributions to Vu between the 
72A-SC(Cd°) and the 72A-SC[Cd-(2e)] is not negligible 
and is caused by the different positions of the Cd nearest 
oxygen neighbors in each case. 



D. Accuracy of the present study: Further 
relaxations and other tests 



The main sources of error of the present study in order 
to compare with experiment are the size of the SC consid- 
ered, the size of the basis, the LDA approximation used 
for the exchange-correlation potential and the relaxation 
process that has been restricted to the six nearest oxy- 
gen neighbors of the Cd atom. To check the accuracy of 
the present study we have performed several additional 
calculations. 



1. Charged impurity: Cd 

We first focus on the case of the charged impurity 
where we have performed the most intensive tests. 

a. Charge of the cell The differences found be- 
tween the results for the 72A-SC[Cd~(2e)] and 72A- 
SC[Cd~ (fluorine)] (Tables HI and HT|) could be taken as a 
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TABLE III: p and d valence contributions to the electric-field gradient at Cd in TiCh, in units of 10 21 V/m 2 , for the neutral 
and charged states of Cd in the 72A-SC. In the last row we give the difference between the values corresponding to Cd~(2e) 
and Cd°. TOT refers to the total (p+d+s-d) valence EFG. 





Vxx 


P 
Vyy 
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Vyy 


Vzz 


Vxx 


TOT 

Vyy 


Vzz 


72A-SC(Cd u ) 

72A-SC[Ccr(2e)] 

diff. 


-1.62 
-2.62 
-1.00 


+3.17 
+3.55 
+0.38 


-1.55 
-0.93 
+0.62 


-6.44 
-0.33 
+6.11 


+4.29 
+ 1.18 
-3.11 


+2.15 
-0.85 
-3.00 


-7.97 
-2.86 
+5.13 


+7.31 
+4.58 
-2.73 


+0.66 
-1.72 
-2.38 



measure of the error performed in simulating the charged 
impurity state. But in fact, we expect the procedure us- 
ing the 72A-SC[Cd~(2c)] to give better results that the 
one using the 72A-SC[Cd~ (fluorine)] since in the former 
the global properties of the system are only smoothly af- 
fected (background density is 0.003 e/A 3 for this SC). 
Moreover, when considering relaxations of atoms more 
distant from Cd than its NN, the presence of fluorine 
atoms will spuriously influence the results, so, in what 
follows we refer to the 72A-SC(-2). 

b. Size effects In order to check how appropriate are 
the dimensions of the 72A-SC (2a x 2b x 3c) used in the 
present work we have performed self-consistent electronic 
structure calculations for a 48A-SC (2a x 2b x 2c) and 
a 96A-SC (2a x 2b x 4c). In these calculations we put 



the nearest oxygen neighbors of Cd at relative positions 
from Cd that correspond to equilibrium in the 72A-SC. 
The size of the forces that oxygen atoms 01 and 02 ex- 
periment in the 48A-SC and 96A-SC is a measure of the 
convergence of the present calculation. We obtained that 
in the 96A-SC the forces on 01 and 02 atoms point out- 
wards (respect to Cd) and are around 0.15 eV/A. Forces 
of this size produce, during the relaxation process in the 
72A-SC, changes in distances smaller than 0.01 A and 
changes in EFG of about 0.2 x 10 21 V/m 2 , so, similar 
changes would be expected if relaxation in the 96A-SC 
was performed. In the case of the 48A-SC forces point 
inwards and are of the same magnitude for 02 atoms 
but they are about 0.65 eV/A in 01 atoms. We see that 
there is a size effect that makes relaxations to be larger 
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for larger SCs, but for the 72A-SC the effect is quite 
small and no significant variations should be expected if 
relaxations in bigger SCs were considered. 

c. Basis size We compute the self-consistent elec- 
tronic structure of the 72A-SC[Cd~(2e)] for the relaxed 
positions of Table Q] increasing the size of the basis to 
7100 LAPW functions (RK M ax=7). Forces on 01 and 
02 atoms are below the tolerance value of 0.025 eV/A 
indicating that the result of the relaxation performed is 
unaltered for a substantial increase of the basis size. In 
fact, the forces on all the other atoms in the cell have 
the same values than using HKmax—Q (within the toler- 
ance) except for the Ti neighbors of Cd in the z direction 
where the forces differ in 0.12 eV/A. Changes obtained 
in Vu components of the EFG tensor are smaller than 
0.05 x 10 21 V/m 2 . We have also considered the inclu- 
sion of a Cd-Ad LO to improve linearization. When a 
Cd-4d LO is introduced with energy at the Fermi level in 
order to improve the description of the impurity states 
no influence in the forces is detected. The change in Vu 
components of the EFG is smaller than 0.1 x 10 21 V/m 2 . 

d. Exchange-correlation potential We performed 
electronic self-consistent calculations and relaxation of 
Cd nearest oxygen neighbors for the system 72A- 
SC[Cd~(2e)] using the generalized gradient approxima- 
tion (GGA) 30 instead of LDA. With the use of this 
parametrization for the exchange-correlation potential 
we obtained 2.18 A and 2.10 A for Cd-Ol and Cd- 
02 distances, respectively (i.e. only a small change of 
0.01 A in Cd-02 distance). For the EFG we obtained 
V 33 =y rF =+4.94 x 10 21 V/m 2 and 77=0.37. 

e. Further relaxations In order to study the effect 
of relaxing the coordinates of atoms beyond the nearest 
neighbors of the Cd atom we have performed the follow- 
ing two additional relaxations: (a) we allow to relax the 
coordinates of all atoms within a cutoff radius i?c=4 A 
centered at Cd (this involves 24 atoms) and (b) idem for 
Rc=4:.6 A (this involves 42 atoms). For radius larger 
than 4.6 A, atoms to be relaxed would be nearer to the 
images of Cd atom from neighboring cells than to the Cd 
itself, so we consider this radius as a limit for the present 
SC. In Table HVl we compare the results obtained for both 
relaxations with the ones corresponding to the NN re- 
laxation. We observe that there is not any qualitative 
change in the results already discuss for the NN relax- 
ation. There are, however, some small variations in the 
values predicted for the EFG. The differences between 
results from relaxations (a) and (b) (that are as larger 
as the differences between relaxations (a) and NN) are 
in part caused because in (b) the relaxation of 03 atom 
is allowed. Atom 03 is the NN of 02 atom in the X 
direction and its relaxation of about 0.04 A allows a fur- 
ther relaxation of about 0.012 A of the 02 atom. The 
fact that atom 03 is the one that experiments the largest 
relaxation (not considering 01 and 02) shows that direc- 
tional bonding plays an important role in this structure. 

In summary, we observed that none of the factors con- 
sidered influence qualitatively the results, but their effect 



TABLE IV: Results of the different relaxations performed for 
the 72A-SC[Cd~ (2e)]. In each case all the coordinates of the 
Na atoms within a radius Rc (in A) are relaxed until forces 
on them are below 0.025 eV/A. All units as in Table I and II. 
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Vzz 


V 


NN 


2.5 


6 
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2.111 


-2.87 


+4.55 


-1.68 


0.26 
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4.0 


24 


2.176 


2.104 


-3.25 


+4.99 


-1.74 


0.30 


(b) 


4.6 


42 


2.187 


2.116 


-3.17 


+4.86 


-1.69 


0.30 



is neither negligible. We therefore confirm our prediction 
about sign (positive) and direction (Y) of V33, but it 
is not possible to perform an exact prediction about its 
magnitude. Ours checks shows that d(Cd-Ol) and d(Cd- 
02) are converged within 0.01 A and V33 and 77 within 
0.5 x 10 21 V/m 2 and 0.1, respectively. Taking the val- 
ues from the last row of Table HVl as our predicted values 
for Vu we obtain discrepancies with experiment ! 16 ' 26 of 
0.37 x 10 21 V/m 2 for | V33 1 and 0.12 for 77 that could be 
attributed to precision errors of the present calculations. 



2. Neutral impurity: CS 

Let us briefly discuss the case of the 72A-SC(Cd°). 
Due to the factor four that exists because of ^-sampling, 
calculations are much more time consuming in this case. 
Hence, we have checked the size of the SC performing 
calculations only on a 48A-SC, and also checked the size 
of the basis repeating the calculation for the 48A-SC with 
KKmax—7- Comparison of the values obtained for EFG 
and forces let us conclude that errors in the 72A-SC(Cd°) 
are expected to be of the same magnitude than the ones 
obtained for the 72A-SC[Cd~(2e)]. Only the effect of 
adding a Cd-4ci LO and the use of GGA instead LDA 
produce a larger variation of the EFG in this case than 
in the charged cell, and is understandable because of the 
larger d contribution to Vu in the case of the neutral 
state of the impurity. Calculations in the 72A-SC(Cd°) 
with this parametrization for the exchange-correlation 
potential and introducing Cd-4d LO predict values of 
V 33 =V X x=-7.90x 10 21 V/m 2 and 77=0.80. In addition, 
we found a force of 0.04 eV/A inwards at 01 atoms, so 
a small refinement in 01 position would be expected. 
About further relaxations, we have only performed relax- 
ation (a) (R c =i A) for the 72A-SC(Cd°) and also obtain 
variations of the same order than in 72A-SC[Cd~(2e)], in 
particular we obtained V33=Vxx=-6.8x 10 21 V/m 2 and 
77=0.97. Then, the same conclusions about accuracy than 
in the 72A-SC[Cd~(2e)] hold for the 72A-SC(Cd°) but 
note that, due to the high value of 77 in this case, the sign 
and direction of eq could change because of precision. 



9 



(a) 



U 



c 

in 

+ 



> 

n 

o 



2.28 

2.22 

2.16-1 

2.10 

2.04 

1.98 

1.92 



-1.75 

-1.80 

-1.85 

-1.90-1 

-1.95 



-2- 
-4 
-6 



d (Cd-02) 



rf(Cd-Ol) 



-n-n-n-n' 



~A- A _ 



A- A 



A 



A A -A- AAA 



ZZ 

0.00 0.03 0.06 0.09 0.12 



(b) 



T3 

u 



> 



O 



2.20 

2.15 

2.10 

2.05 

2.00 

1.95 
6 

5 

4 

3 

2 

1 



6 
4 
2 

-2 
-4 
-6 



rf(Cd-Ol) 
• • • • • 



d (Cd-02) 



F(01) 




F(02) 



• * . 



zz- 



XX 



Ad./d. 



0123456789 
Step of relaxation 



FIG. 8: (a) Isotropic relaxation process. Distances di=d(Cd-Oi) and a!2=d(Cd-02), total energy, and Vu as a function of the 
relative displacement, Adi/di, of the nearest O atoms for the 72A-SC[Cd~ (2e)] system (Adi — di — di,d® is di for the pure 
system); (b) Free relaxation process. Distances d\ and efe, forces at Ol and 02, and Vu as function of the step of the relaxation 
in the 72A-SC[Cd-(2e)] system. 



E. Comparison with other calculations 

The simplest and most widely used approximation for 
the calculation of the EFG at a probe-atom is the point- 
charge model (PCM). 2 In this approximation the EFG 
tensor at the probe site is (1 — 7oo)V^ att ', where V^ u - 
is the EFG tensor produced by valence nominal charges 
located at the ion positions in the lattice, and 7oo is the 
Sternheimer antishiclding factor— that depends only on 
the probe-atom. In this way, this model assumes that 
the symmetry and orientation of the EFG tensor at im- 
purity sites are unaltered by its presence. PCM gives 
forTi0 2 (Cd): V 33 =V XX =-2.27 x 10 21 V/m 2 and n= 0.40 
when a value of -29.27 is used for 700^ In Ref. [l8| the 
authors assume that relaxation is responsible for the dis- 
agreement of PCM predictions with experiment and spec- 
ulate that an isotropic relaxation of 0.04 A outwards 
of all the Cd nearest oxygen neighbors would produce 
the desired result for | V33 1 and i]. Our results indi- 
cate that relaxations are not isotropic and they are so 



much larger than this, but even if the relaxed coordi- 
nates from our calculation (see Table I) were used, the 
PCM would fail in the description of the EFG giving 
V33=Vxx= -8.20 x 10 21 V/m 2 in clear contradiction with 
V 33 =V Y y= +4.55 x 10 21 V/m 2 that we obtain from the 
self-consistent FLAPW calculation. The disagreement of 
PCM with the FLAPW prediction could not either be at- 
tributed in this case to a change in the value of 700 since 
the sign and directions of both predictions are different. 
Thus, it is clear that the problem of EFG at cationic sites 
in Ti0 2 is too complicated to be described even approx- 
imately by simple PCM calculations. 

In our previous work 2 ^ we performed self-consistent 
electronic structure calculations of this system with a 
twelve-atom super-cell (12A-SC: a x b x 2c). Relax- 
ations were performed only for the neutral charge state 
of the impurity and at the end of the relaxation pro- 
cess two electrons were added to compute EFG without 
computing the self-consistent potential of the charged 
cell. The structural relaxation obtained by this proce- 
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dure were smaller than in the present work but account 
for the inversion of Cd-Ol and Cd-02 distances [d(Cd- 

01) =2.12 A, d(Cd-O2)=2.07 A] with respect to the un- 
relaxed structure (see Table IJ). The description obtained 
for the EFG was very similar to the one of the 72A- 
SC[Cd~(2e)], but we understand that this agreement is 
somehow fortuitous because usage of 12A-SC relaxed co- 
ordinates to compute the EFG in the 72A-SC[Cd"(2e)] 
gives very different results. 

In Ref. [Til the authors have used a very similar ap- 
proach as that we used in Ref. using a 12A-SC, 
but they assumed that the relaxations of the nearest 
oxygen neighbors of the Cd atom were isotropic. As 
a consequence, they obtained d(Cd-Ol)=2.04 A, c?(Cd- 

02) =2.08 A, and also a very different result for EFG 
(see Table I), V ss =Vzz=-5.Q9 x 10 21 V/m 2 and ?7=0.39 
(for a carefully comparison of these two calculations see 
Ref HI). 

A question that arises at this point is if self-consistent 
electronic FLAPW calculations performed with a con- 
verged SC give EFG compatible with experience or not 
in the case that relaxations of the oxygen NN of the Cd 
atom are constrained to be isotropic. This is an inter- 
esting point in order to know if available experimental 
data is enough to refute the assumption that the oxygen 
NN relax isotropically. We have performed self-consistent 
calculations for the system 72A-SC[Cd~(2e)] for differ- 
ent positions of 01 and 02 atoms but moving them 
outwards, keeping the relation d(Cd-01)/d(Cd-02) con- 
stant. We determined the equilibrium position of oxy- 
gen atoms as the one that produce a minimum in the 
energy [see Fig. Ufa)]. We obtain that Cd-0 distances 
are 2.12 A and 2.16 A for 01 and 02 atoms, respec- 
tively [a relaxation of 9% of the unrelaxed distances, see 
Fig- EI a)]. If we compare Figs. [5{a) and (b), we can note 
that due to the assumption of isotropic relaxations, there 
is no inversion in the Cd-0 distances and, in consequence, 
the strong change in the EFG components does not take 
place. In particular, there is not change in sign and ori- 
entation of V33 as in the case of our free relaxation. At 
the equilibrium position we obtained for V33 a value of - 
4.46 x 10 21 V/m 2 (pointing in [001]-direction) and a high 
77 value of 0.91, confirming that an isotropic relaxation is 
not consistent with the experimental data. 

IV. CONCLUSIONS 

In this work we have studied through a series of 
first-principles calculations the problem of a Cd impu- 
rity substitutionally-located at the cationic site in ru- 
tile Ti02- The main result of our work, i.e. that Cd 
introduces in the host fairly anisotropic relaxations of 
its nearest oxygen neighbors and that this produces a 
change of orientation of V33 from the [001] to the [110] 
direction when a Ti atom is replaced by a Cd atom in 
pure Ti02, was briefly presented in a recent work with 
the experimental confirmation of the last prediction 16 . 



In this work we have presented details about the elec- 
tronic structure of the different impurity systems consid- 
ered. We have considered atomic relaxations and elec- 
tronic structure self-consistently and have obtained that 
both aspects of the problem interact with each other. 
We obtained that atomic relaxations are different for the 
charged and neutral state of the impurity and that, on 
the other side, the relaxation process produces a dras- 
tic variation in the asymmetry of the charge distribution 
near the probe-atom for a given charge state of the im- 
purity, which is detected in the strong variation of the 
EFG tensor. We have shown that the huge difference in 
the values of the asymmetry parameter 77 between the 
charged and neutral state of the impurity arises because 
of the filling of the impurity level at the Fermi energy. 
This difference in the rj value determines, through com- 
parison with experiment, that Cd is in a charged state 
when it is introduced as impurity in Ti02 at room tem- 
perature. From these results we have confirmed that the 
EFG tensor is a very useful magnitude because it is sen- 
sible to subtle details of the electronic structure and it 
can be determined experimentally with high resolution. 
We have performed a series of checks of the accuracy of 
the present calculations in order to show that all the pre- 
dictions of this work are the same if an increment in the 
basis size, k-mesh or size of the SC are considered or if a 
different exchange-correlation potential is used. We have 
also shown that considering relaxations beyond nearest 
neighbors does not produce any qualitative change in our 
results. Finally, we have checked that the hypothesis of 
isotropic relaxations and the use of PCM approximations 
give results incompatible with experiments and with our 
calculations. From our results it is clear that the problem 
of the EFG at Cd impurities in Ti02 is too complicated 
to be described (even approximately) by simple models 
like PCM, antishielding factors and isotropic relaxations. 
We can conclude that a proper theoretical description of 
electronic properties of metal impurities in oxide semi- 
conductors should consider self-consistently the charge- 
state of the impurity and the impurity-induced distor- 
tions in the host, specially in the first shell of neighbors 
of the impurity. 
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